Inelastic collisions in molecular nitrogen at low temperature (2<T<50 K)

Theory and experiment are combined in a novel approach aimed at establishing a set of two-body state-to-state rates for elementary processes ij->lm in low temperature N2:N2 collisions involving the rotational states i, j, l, m. First, a set of 148 collision cross sections is calculated as a function of the collision energy at the converged close-coupled level via the MOLSCAT code, using a recent potential energy surface for N2–N2. Then, the corresponding rates for the range of 2<T<50 K are derived from the cross sections. The link between theory and experiment, aimed at assessing the calculated rates, is a master equation which accounts for the time evolution of rotational populations in a reference volume of gas in terms of the collision rates. In the experiment, the evolution of rotational populations is measured by Raman spectroscopy in a tiny reference volume 2E-3 mm3 of N2 traveling along the axis of a supersonic jet. The calculated collisional rates are assessed experimentally in the range of 4<T<35 K by means of the master equation, and then are scaled by averaging over a large set of experimental data. The scaled rates account accurately for the evolution of the rotational populations measured in a wide range of conditions. Accuracy of 10% is estimated for the main scaled rates.This work has been supported by the Spanish Ministerio de Educación y Ciencia, research Project Nos. FIS2004-02576, HF2004-232, ESP2004-21060-E, and ASTROCAM network. J.P.F. is indebted to the CSIC for an I3P grant.Peer reviewe

The Abundance of SiC2 in Carbon Star Envelopes: Evidence that SiC2 is a gas-phase precursor of SiC dust

Silicon carbide dust is ubiquitous in circumstellar envelopes around C-rich AGB stars. However, the main gas-phase precursors leading to the formation of SiC dust have not yet been identified. The most obvious candidates among the molecules containing an Si--C bond detected in C-rich AGB stars are SiC2, SiC, and Si2C. We aim to study how widespread and abundant SiC2, SiC, and Si2C are in envelopes around C-rich AGB stars and whether or not these species play an active role as gas-phase precursors of silicon carbide dust in the ejecta of carbon stars. We carried out sensitive observations with the IRAM 30m telescope of a sample of 25 C-rich AGB stars to search for emission lines of SiC2, SiC, and Si2C in the 2 mm band. We performed non-LTE excitation and radiative transfer calculations based on the LVG method to model the observed lines of SiC2 and to derive SiC2 fractional abundances in the observed envelopes. We detect SiC2 in most of the sources, SiC in about half of them, and do not detect Si2C in any source, at the exception of IRC +10216. Most of these detections are reported for the first time in this work. We find a positive correlation between the SiC and SiC2 line emission, which suggests that both species are chemically linked, the SiC radical probably being the photodissociation product of SiC2 in the external layer of the envelope. We find a clear trend in which the denser the envelope, the less abundant SiC2 is. The observed trend is interpreted as an evidence of efficient incorporation of SiC2 onto dust grains, a process which is favored at high densities owing to the higher rate at which collisions between particles take place. The observed behavior of a decline in the SiC2 abundance with increasing density strongly suggests that SiC2 is an important gas-phase precursor of SiC dust in envelopes around carbon stars.Comment: Published in A&A. 16 pages and 10 figure

Clues to NaCN formation

ALMA is providing us essential information on where certain molecules form. Observing where these molecules emission arises from, the physical conditions of the gas, and how this relates with the presence of other species allows us to understand the formation of many species, and to significantly improve our knowledge of the chemistry that occurs in the space. We studied the molecular distribution of NaCN around IRC +10216, a molecule detected previously, but whose origin is not clear. High angular resolution maps allow us to model the abundance distribution of this molecule and check suggested formation paths. We modeled the emission of NaCN assuming local thermal equilibrium (LTE) conditions. These profiles were fitted to azimuthal averaged intensity profiles to obtain an abundance distribution of NaCN. We found that the presence of NaCN seems compatible with the presence of CN, probably as a result of the photodissociation of HCN, in the inner layers of the ejecta of IRC +10216. However, similar as for CH 3 CN, current photochemical models fail to reproduce this CN reservoir. We also found that the abundance peak of NaCN appears at a radius of 3 x 10 15 cm, approximately where the abundance of NaCl, suggested to be the parent species, starts to decay. However, the abundance ratio shows that the NaCl abundance is lower than that obtained for NaCN. We expect that the LTE assumption might result in NaCN abundances higher than the real ones. Updated photochemical models, collisional rates, and reaction rates are essential to determine the possible paths of the NaCN formation.Comment: 7 pages, 10 figures. Accepted for publication in A&A letter

The abundances of polyacetylenes towards CRL618

We present a mid-infrared high spectral resolution spectrum of CRL618 in the frequency ranges 778-784 and 1227-1249 cm^-1 (8.01-8.15 and 12.75-12.85 um) taken with the Texas Echelon-cross-Echelle Spectrograph (TEXES) and the Infrared Telescope Facility (IRTF). We have identified more than 170 ro-vibrational lines arising from C2H2, HCN, C4H2, and C6H2. We have found no unmistakable trace of C8H2. The line profiles display a complex structure suggesting the presence of polyacetylenes in several components of the circumstellar envelope (CSE). We derive total column densities of 2.5 10^17, 3.1 10^17, 2.1 10^17, 9.3 10^16 cm^-2, and < 5 10^16 cm^-2 for HCN, C2H2, C4H2, C6H2, and C8H2, respectively. The observations indicate that both the rotational and vibrational temperatures in the innermost CSE depend on the molecule, varying from 100 to 350 K for the rotational temperatures and 100 to 500 K for the vibrational temperatures. Our results support a chemistry in the innermost CSE based on radical-neutral reactions triggered by the intense UV radiation field.Comment: 9 pages, 4 figures, 1 table; accepted for publication in The Astrophysical Journa

High-J v=0 SiS Maser Emission in IRC+10216: A New Case of Infrared Overlaps

We report on the first detection of maser emission in the J=11-10, J=14-13 and J=15-14 transitions of the v=0 vibrational state of SiS toward the C-rich star IRC+10216. These masers seem to be produced in the very inhomogeneous region between the star and the inner dust formation zone, placed at 5-7 R*, with expansion velocities below 10 km/s. We interpret the pumping mechanism as due to overlaps between v=1-0 ro-vibrational lines of SiS and mid-IR lines of C2H2, HCN and their 13C isotopologues. The large number of overlaps found suggests the existence of strong masers for high-J v=0 and v=1 SiS transitions, located in the submillimeter range. In addition, it could be possible to find several rotational lines of the SiS isotopologues displaying maser emission.Comment: 4 pages, 1 figure, published in the ApJ Letter

Detection of the S(1) Rotational Line of H2 toward IRC+10216: A Simultaneous Measurement of Mass-Loss Rate and CO Abundance

We report the first detection of the S(1) pure rotational line of ortho-H2 at 17.04 um in an asymptotic giant branch star, using observations of IRC+10216 with the Echelon-cross-echelle Spectrograph (EXES) mounted on the Stratospheric Observatory for Infrared Astronomy (SOFIA). This line, which was observed in a very high sensitivity spectrum (RMS noise ~0.04% of the continuum), was detected in the wing of a strong telluric line and displayed a P Cygni profile. The spectral ranges around the frequencies of the S(5) and S(7) ortho-H2 transitions were observed as well but no feature was detected in spectra with sensitivities of 0.12% and 0.09% regarding the continuum emission, respectively. We used a radiation transfer code to model these three lines and derived a mass-loss rate of 2.43(0.21)E-05 M_sun/yr without using the CO abundance. The comparison of this rate with previous estimates derived from CO observations suggests that the CO abundance relative to H2 is 6.7(1.4)E-04. From this quantity and previously reported molecular abundances, we estimate the O/H and C/H ratios to be 3.3(0.7)E-04 and >5.2(0.9)E-04, respectively. The C/O ratio is >1.5(0.4). The absence of the S(5) and S(7) lines of ortho-H2 in our observations can be explained by the opacity of hot dust within 5R* from the center of the star. We estimate the intensity of the S(0) and S(2) lines of para-H2 to be ~0.1% and 0.2% of the continuum, respectively, which are below the detection limit of EXES.Comment: 7 pages, 4 figures, accepted for publication in ApJ Letter

Multi-frequency high spectral resolution observations of HCN toward the circumstellar envelope of Y CVn

High spectral resolution observations toward the low mass-loss rate C-rich, J-type AGB star Y CVn have been carried out at 7.5, 13.1 and 14.0 um with SOFIA/EXES and IRTF/TEXES. Around 130 HCN and H13CN lines of bands v2, 2v2, 2v2-v2, 3v2-2v2, 3v2-v2, and 4v2-2v2 have been identified involving lower levels with energies up to ~3900 K. These lines have been complemented with the pure rotational lines J=1-0 and 3-2 of the vibrational states up to 2v2 acquired with the IRAM 30 m telescope, and with the continuum taken with ISO. We have analyzed the data with a ro-vibrational diagram and a code which models the absorption and emission of the circumstellar envelope of an AGB star. The continuum is produced by the star with a small contribution from dust grains comprising warm to hot SiC and cold amorphous carbon. The HCN abundance distribution seems to be anisotropic. The ejected gas is accelerated up to the terminal velocity (~8 km/s) from the photosphere to ~3R* but there is evidence of higher velocities (>9-10 km/s) beyond this region. In the vicinity of Y CVn, the line widths are as high as ~10 km/s, which implies a maximum turbulent velocity of 6 km/s or the existence of other physical mechanisms probably related to matter ejection that involve higher gas expansion velocities than expected. HCN is rotationally and vibrationally out of LTE throughout the whole envelope. A difference of about 1500 K in the rotational temperature at the photosphere is needed to explain the observations at 7.5 and 13-14 um. Our analysis finds a total HCN column density that ranges from ~2.1E+18 to 3.5E+18 cm^{-2}, an abundance with respect to H2 of 3.5E-05 to 1.3E-04, and a 12C/13C isotopic ratio of ~2.5 throughout the whole envelope.Comment: 24 pages, 11 figures, 3 tables, accepted for publication in A&